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Abstract. Magnetic confinement in graphene has been of recent and growing interest because its 
potential applications in nanotechnology. In particular, the observation of the so called magnetic 
edge states in graphene has opened the possibility to deepen into the generation of spin currents 
and its applications in spintronics. We study the magnetic edge states of quasi-particles arising 
in graphene monolayers due to an inhomogeneous magnetic field of a magnetic barrier in the 
formalism of the two-dimensional massless Dirac equation. We also show how the solutions of such 
states in each of both triangular sublattices of the graphene are related through a supersymmetric 
transformation in the quantum mechanical sense. 
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INTRODUCTION 

Graphene is a novel two dimensional material that has opened a new bridge of common 
interests between the condensed matter and high energy physics communities. It con- 
sists of a single atomic monolayer of carbon atoms arranged in a honeycomb hexagonal 
crystal lattice that can be described through two triangular Bravias sublattices A and B. 
Some of its properties of interest for the nanoelectronics are its high electrical and ther- 
mal conductivity as well as its high elasticity and toughness, making it a lightweight and 
flexible material that could be used instead of silicon in some semiconductor electronic 
devices. The structure of the energy bands of graphene is such that the dispersion relation 
is a linear one, implying that the charge carriers have a zero effective mass. Then, the 
graphene can be studied through the Dirac equation in (2 + 1)D in the zero mass limit, 
where the two irreducible representations of the Dirac matrices describe two species of 
particles, one for each triangular sublattice and labeled by a pseudo-spin originated by 
their interactions with the lattice [1]. Interacting with an external magnetic field, the 
electrons in the crystal lattice of graphene are subjected to an effective potential that has 
a minimum below the zero of the external potential [2]. This implies the existence of 
bound states solutions, giving rise to some phenomena of confinement like the so called 
magnetic edge states, magnetic quantum dots and magnetic spin currents, all of them 
of growing interest in the nanotechnology. Here we obtain the magnetic edge states in 
graphene for a static external magnetic barrier making use of the supersymmetric char- 
acter, in the quantum mechanical sense (SUSY-QM), of the Dirac equation for external 
and time-independent magnetic fields [3]. 



DIRAC EQUATION IN (2+l)D AND SUSY-QM 



In the context of the relativistic quantum mechanics, a super symmetric Hamiltonian 
can be written in the general form as //susy = {Q, Q?}, with the supercharges Q and 
Q} being nilpotent operators. We are interested in the (2 + 1)D Dirac equation in the 
massless limit, 

(y.n)¥(v) = 0, (1) 

with = id^ +eAf 1 l and *P(x,f) = \i/(x)e~ lEt . In this case, there are two irreducible 
representations for the 7 matrices given in terms of the 2x2 Pauli matrices. We will work 
with the Jackiw representation, for which the first irreducible representation, labeled 
as A, corresponds to the set {y° = 03, y 1 = iO\,y 2 = 1O2}, and the second irreducible 
representation, B, is given by {y° = 03, y 1 = io\ , y 2 = — io{\. 

Taking the square of the (y-IT) operator in Eq. (1), (y-TI) 2 ^, a time-independent 
Schrodinger-type equation for i/a(x) arises in terms of a SUSY-QM Hamiltonian, H a R = 

{Qr-> Or} = — (7 - n) 2 + riQ, which depends on both, the spin of the particles (a) and 
on the irreducible representation (R = {A, B}) of the Dirac matrices, as H aj R\j/(x) — 
£ ,2 i//(x). The supercharges, in each irreducible representation, are Qr = (Jl\ + U2)Or 

and Q\ = (111 — n2)<7^", with = (c 1 ±io 2 )/2 = o^. There are two corresponding 
effective potentials V^, one per each irreducible representation, called SUSY partner 
potentials. For an external static magnetic field described in the Landau gauge by the 
tri-potential A^ = (0,0,W(x)), these are given in terms of the superpotential W = 
eW(x) + as V a ,R = W + ( — l) R oW' with W = dW(x) jdx, P2 being the component 
of the momentum along the y direction and R = 1 , 2 for the irreducible representation A 
and B respectively. 

As a consequence of the SUSY-QM character of H aR , the Dirac equation (1) for 
any external static magnetic field, reduces to a Pauli one with effective mass m = 
1/2 and gyromagnetic ratio g = 2 [5]. Due to the fact that there are two irreducible 
representations for the Dirac matrices in the (2 + 1)D case, there is also a direct relation 
between the solutions of the wave function for different spin eigenvalues, which shows 
itself through the operators. 

MAGNETIC BARRIER 

Applying the formalism of SUSY-QM for a magnetic barrier of length d perpendicular 
to the plane of a graphene monolayer, given by B(x) = Bo&(d 2 —x 2 ), each one of the 
effective potentials V g _r corresponds to a truncated harmonic oscillator that in general is 
not symmetric. The condition for the existence of bound states solutions, referred in the 
literature as edge states (because they are localized at the borders of the external applied 
potential), in this case is given if the energy E of the charge carries satisfies the relation 



We work in natural units such that % = c = 1 , being c the Fermi velocity, which for graphene plays the 
role of the speed of light and is two orders of magnitude smaller than c (c ~ c/300) [4]. 



£ 2 B E 2 < V ( s \ being s = sign(/> 2 40, V (±) = (p 2 £B±d/£ B ) and 4 = ^Jhc/{eB) is the 
magnetic length. 

Figure 1 shows the probability density function | ^(x/ig) | 2 for the calculated magnetic 
edge states for a magnetic barrier of length d = 3£b for two values of the £ B p 2 parameter. 
The spin direction has been selected as(— 1) R <J = — 1. 




FIGURE 1. The probability density function \\i/(x/£b)\ 2 of the magnetic edge states for a magnetic 
barrier of length d = 3£b for (— l) R o = —1 and pits =-1.0 (left), 0.0 (right). The magnetic barrier, the 
effective truncated harmonic oscillator, as well as a zoom of this last in the energy region of the bound 
states, are schematically shown in both cases. 



CONCLUDING REMARKS 

SUSY-QM applied to graphene in a static and inhomogeneous external magnetic field, 
allows to map the wave functions of the charge carriers in one of the two triangular 
sublattices, with definite induced pseudo-spin, into the other. In this way the solutions 
of the (2 + 1)D Dirac equation in each of the two irreducible representations are closely 
related. In the zero mass limit, the relativistic description of the system has been reduced 
to a non-relativistic one characterized through a time-independent Schrodinger-type 
equation. Under this formalism, the magnetic edge states of graphene for a magnetic 
barrier were obtained. 



REFERENCES 

1. G. W. Semenoff, Phys. Rev. Lett. 53, 2499-2452 (1984); D. P. DiVicenzo, and E. J. Mele, Phys. Rev. 
B29, 1685-1694(1984). 

2. A. De Martino, L. Dell' Anna, and R. Egger, Phys. Rev. Lett. 98, 066802-1-4 (2007); Solid State 
Commun. 144, 547-550 (2007); M. Ramezani Masir, A. Matulis, and F. M. Peeters, Phys. Rev. B 79, 
155451-1-8(2009). 

3. F. Cooper, A. Khare, and U. Sukhatme, Supersymmetry in Quantum Mechanics, World Scientific, 
Singapore, 2001. 

4. A. Geim, and K. Novoselov, Nature Mater. 6, 183-191 (2007); N. M. R. Peres, F. Guinea, and A. H. 
Castro Neto, Phys. Rev. B 73, 12541 1-1-23 (2006). 

5. M. Moreno, and R. M. Mendez-Moreno, "Generalized supersymmetric quantum mechanics," in Pro- 
ceedings of the Workshop on High-energy Phenomenology, edited by M. Perez and R. Huerta, World 
Scientific, Singapore, 1992, pp. 365-375. 



